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We demonstrate enhanced di-Higgs production at the LHC in the presence of modifications of 
the effective couplings of Higgs to gluons from new, light, colored scalars. While our results apply 
to an arbitrary set of colored scalars, we illustrate the effects with a real color octet scalar - a 
simple, experimentally viable model involving a light (~ 125-300 GeV) colored scalar. Given the 
recent LHC results, we consider two distinct scenarios: First, if the Higgs is indeed near 125 GeV, 
we show that the di-Higgs cross section could be up to nearly 10 times the Standard Model rate 
for particular octet couplings and masses. This is potentially observable in single Higgs production 
modes, such as pp — > hh — > jjbb as well as pp — > hh — > r + T~bb where a small fraction of the 77 
or t + t~ events near the putative Higgs invariant mass peak contain also a bb resonance consistent 
with the Higgs mass. Second, if the Higgs is not at 125 GeV (and what the LHC has observed 
is an impostor), we show that the same parameter region where singly-produced Higgs production 
can be suppressed below current LHC limits, for a heavier Higgs mass, also simultaneously predicts 
substantially enhanced di-Higgs production. We point out several characteristic signals of di-Higgs 
production with a heavier Higgs boson, such as pp — > hh — > W + W~W + W~ , which could use 
same-sign dileptons or trileptons plus missing energy to uncover evidence. 



I. INTRODUCTION 

Di-Higgs production in the Standard Model (SM) has 
a very small rate at the LHC [DH]. For mh = 125 GeV, 
the leading order cross section is a(pp — > hh) = 4 (16) fb 
at y/s = 8 (14) TeV. This is much smaller than the single 
Higgs production cross section due to the larger partonic 
energy needed to produce two Higgs bosons, as well as 
an accidental cancellation between the s-channel and box 
diagram contributions to the amplitude gg — > hh. Hence, 
like other accidentally suppressed processes in the Stan- 
dard Model, di-Higgs production provides a great oppor- 
tunity for new physics to be observed. In this paper we 
demonstrate that in the presence of a a general set of 
light colored scalars, the di-Higgs production cross sec- 
tion can be enhanced by orders of magnitude above the 
Standard Model rate. Light colored scalars that couple to 
the Higgs boson are well known to have a profound effect 
on single Higgs production [5"Ml3j. Di-Higgs production 
has been considered previously in several different con- 
texts [I4H36] . 

Since direct searches at the LHC for light colored 
scalars that couple to the Higgs have many constraints, 
depending on the model and the decay modes, the bur- 
den is on us to provide a "benchmark model" that evades 
these constraints while providing the enhancements in di- 
Higgs production that we find. Several possible represen- 
tations could be considered, including the superpartners 
to the top, vector-like sets of particles (that may or may 
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not carry electroweak quantum numbers), or real repre- 
sentations such as an electroweak neutral color octet. 

A real scalar octet is an interesting example of a col- 
ored scalar that can be effectively hidden in the LHC 
data. In the model that we consider, these scalars are 
pair-produced (up to small corrections coming from loop 
processes), and decay through loop-level processes to a 
pair of gluons. Thus, the signature is two pairs of jets 
with equal invariant mass, making four total jets. AT- 
LAS [37] and CMS [38] have multi-jet searches that are 
potentially sensitive to this scalar; the ATLAS study re- 
stricts M s > 125 GeV, while the CMS study [3B] is sen- 
sitive only for scalars exceeding 320 GeV. This leaves a 
wide range of real scalar color octets are allowed by LHC 
constraints. A real scalar color octet S a can couple to the 
Higgs boson through the renormalizable Higgs-portal in- 
teraction of the form (k/2)S^H' H . Vital to our study is 
the possibility that n can be negative, as we will see. 

The correlation between effects in single-Higgs and di- 
Higgs production has been studied previously in terms 
of contact interactions [26j [29]. For very heavy col- 
ored scalars - meaning scalars whose mass is large com- 
pared to the electroweak scale, Mi ^> v, as well as to 
the characteristic energy of di-Higgs production (roughly 
Mj ^> s ~ O(10 ml) - the effective operators provide 
a reasonable description. As we will see, for the col- 
ored scalars we consider in this paper, between about 
100 to 300 GeV, the full momentum-dependence of the 
one-loop calculation is essential to accurately estimate 
the di-Higgs enhancement. We show this explicitly in 
Appendix [B] 

Finally, it may seem somewhat quixotic that we con- 
sider mh = 125 GeV as well as m/, > 125 GeV, given the 



strong evidence from the LHC for a new particle with 
properties consistent with a Higgs near rrih = 125 GeV 
[3"9l |4T)] . While much work remains to be done to verify 
that the 125 GeV particle is indeed the (or a) Higgs bo- 
son, there is a more direct reason for our continued inter- 
est in a heavier Higgs boson. As we showed in Ref. [12], 
colored scalars can suppress single Higgs production well 
below the current bounds from LHC for a wide range of 
Higgs masses. So the argument that "we have already 
searched for a heavier Higgs boson and did not find it, 
so the Higgs must be the 125 GeV Higgs-like particle" 
is simply wrong. This reasoning is wrong because the 
argument applies only to the Standard Model. The gen- 
eral class of models we consider in this paper - colored 
scalars with Higgs portal couplings - provide a clear class 
of counterexamples. Of course we are not disputing the 
strong evidence for a 125 GeV particle; instead, we be- 
lieve maintaining a healthy dose of skepticism regarding 
the true identity of this particle, given the wide num- 
ber of impostors [H] that could be masquerading as a 
Higgs-like resonance. 



of 1/Mi. The leading interactions of the Higgs to glu- 
ons can be determined by matching the strong coupling 
constant in the low energy and high energy theory [15] , 
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where we have included the top quark and an arbi- 
trary set of colored scalars with explicit Higgs field- 
dependence. For a real scalar field, C, is the Dynkin 
index, e.g. Cj = 3 for a color octet; for a complex repre- 
sentation, replace Cj — > 2Ci. This leads to the one- loop 
effective Lagrangian, 
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II. EFFECTIVE COUPLINGS 

Low energy theorems for Higgs physics provide pow- 
erful methods to determine the effective Higgs couplings 
to both Standard Model particles as well as new physics 
15, 42-46J. We are interested in extending the Stan- 
dard Model to include a general set of colored scalars 
in arbitrary representations of QCD. The multiplicity of 
scalars (number of flavors) is taken into account, while 
our predictions for the single and di-Higgs rates are oth- 
erwise independent of the electroweak quantum numbers 
at leading order in the couplings. 

The minimal Lagrangian for colored scalars in complex 
representations is 

C c = (D^iD^i) - mfyUi - K^l&rfH (1) 



while for real scalars the Lagrangian is 
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We have not included quartic (or possibly cubic) self- 
interactions, nor interactions among different flavors of 
scalars, since these couplings will not play any role 
in our calculations of di-Higgs production. The field- 
independent scalar mass-squared is 
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while the Higgs field-dependent mass is 
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both of which apply to scalars in real or complex repre- 
sentations. When TO/j <C M,, the colored scalars can be 
integrated out, resulting in an effective theory in powers 



where we have shifted Gt v back to the canonical basis. 

Expanding Eq. (JTjj) to first order in h/v, we obtain the 
single Higgs effective interaction with colored scalars, 
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The single Higgs effective interaction, Eq. Q, is for ex- 
ample identical to the result obtained from our previous 
study of Higgs suppression through colored scalars in the 
limit m h < M t [Eq. (2.7) from Ref. [L2]]. 

The 4-point di-Higgs effective interaction is similarly 
obtained, 
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The top quark contributions to the effective couplings 
can also be obtained by expanding Eq. ([6]), 
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There are two critical observations we can make at the 
effective interaction level. First, when Mj > x/|fCj]u, the 
relative sign of the the top quark contribution to the ef- 
fective operator h 2 G a ^ v G a ^ v is is opposite (the same as) 
that of the scalar contribution when m > (/tj < 0). 
Second, in the limit mh,v <C M i; we see that the same 
coefficient that determines the single Higgs effective in- 
teraction also uniquely determines the di-Higgs effective 
interaction. This will be important for gaining some 
qualitative understanding of our di-Higgs cross section 
results. 
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FIG. 1. Effective operator contributions to gg —¥ hh. 



III. CONTRIBUTIONS TO DI-HIGGS 
PRODUCTION 

We now carry out the calculation of di-Higgs produc- 
tion at a hadron collider. Light quarks play no role in 
the leading-order di-Higgs production cross section, so we 
are left with the gluon-induced partonic process gg — > hh. 
We calculate the amplitudes for di-Higgs production in 
two independent mcthonds. First we calculate the am- 
plitudes using the effective interactions derived in Sec.[H] 
This allows us to obtain simple, analytic formulae that 
provide several qualitative features of the full calculation. 
Next, we carry out a full momentum-dependent one- loop 
calculation of gg — > hh including top quarks and an ar- 
bitrary set of colored scalars. 

The amplitude for gg — > hh can be decomposed into 
two non- interfering Lorentz structures following Ref. pQ, 

M(g a gb -> hh)nu = VP^ ah + QQ^S ab , (10) 
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In these formula, p\ and pi are the incoming momenta 
of the gluons, k\ is the outgoing momenta of one of the 
Higgs bosons, and px is the transverse momenta of the 
Higgs: p\ — (ut— m 4 H )/s. The is the only possible 
alternative Lorentz structure once the Ward identities 
and orthogonality to P^ have been imposed. 



A. Di-Higgs Amplitude from Effective Couplings 
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P^y structure. 



The effective couplings given by Eqs. ([7]), ([8]), 

tlead only to contributions to the 
e Feynman diagrams include both the 4-point ef- 
fective interaction h 2 G a ^ v G a4il ' as well as the diagram 
with s-channel Higgs exchange involving the 3-point ef- 
fective interaction hG a ilv G a ' llv and the triple-Higgs self- 
interaction, shown in Fig. [I] We obtain 



where we have neglected the width of Higgs in the s- 
channel propagators. The top quark contributions to the 
amplitude interfere destructively. This, combined with 
the necessity to sample parton distribution functions at 
larger x (to obtain larger s), is the central reason that 
the di-Higgs production cross section at a hadron collider 
is so small in the Standard Model. 

The amplitudes derived from the effective interactions 
allow us to make several interesting qualitative observa- 
tions about the scalar contributions: 

• For Ki > with Mi <; i^/jTcTj?;, the top contribution 
and the scalar contributions destructively interfere, 
weakening the effects of colored scalars on di-Higgs 
production. 

• For ^ < with Mi > y^KTjv, the top contribution 
and the scalar contributions constructively inter- 
fere, strengthening the effects of colored scalars on 
di-Higgs production. 



When Mi < 



Ki\v, the third contribution in 
Eq. ( 14 ) begins to affect the scalar amplitude for di- 
Higgs production. However, momentum-dependent 
corrections proportional to s/M 2 are at least as im- 
portant, since s > irri^. In this region, the 1/Mi 
expansion is no longer valid, and we need the full 
momentum-dependent loop functions to make ac- 
curate quantitative calculations. 

The two main results we find from the di-Higgs ampli- 
tudes calculated with the effective operators are that: i.) 
we expect a considerably larger cross section when Ki < 
0, and Mi J> y^ftTju; ii.) we expect a strong correlation 
between single Higgs production and di-Higgs production 
once Mi is large enough for the effective operators to re- 
produce the full momentum-dependent one- loop results. 



B. Di-Higgs Amplitude at One-Loop 

We now turn to the full one-loop leading-order cal- 
culation of di-Higgs production including the top quark 
and the scalars. The Standard Model contribution to 
gg — > hh comes from top quark triangle diagrams stitched 
to a triple-Higgs vertex via a Higgs boson propagator, as 
well as box diagrams with two top Yukawa coupling in- 
sertions. The full momentum-dependent expressions for 
V, Q for the top contributions can be found in Ref. [T]. 

The scalar octet loops can be similarly classified into 
triangles and boxes (or by powers of the Higgs-portal cou- 
pling k). The complete set of (leading order) diagrams 
are shown below in Fig. [2j The scalar triangle diagram 
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FIG. 2. Scalar loop contributions to gg — > hh. The three diagrams on left are 0{a s n), while the four diagrams on the right 
are 0(a s K 2 ), where k is the Higgs-portal coupling. 



contributions only have P^ v gauge structure: 
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The first term in the parenthesis comes from attaching 
the 4-point vertex in the k interaction to the scalar loop, 



while the second term comes from connecting a Higgs 
propagator and 3-point vertex to the triple-Higgs self- 
interaction. 

The box diagrams involving scalars (as well as top 
quarks) contribute to both (P^, Q^v) Lorentz struc- 
tures. We evaluate the scalar contribution to V and Q 



(as well as "Ptn in Eq. (15)) in terms of the Passarino- 
Vcltman one-loop functions given in Appendix [Aj We 
obtain: 
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Adding the scalar loop V,Q to the top loop contribu- 
tions gives us the total M(gg —> hh)u V amplitude. Squar- 
ing and adding phase space, color- and spin-averaging 
factors, we arrive at the differential partonic cross sec- 
tion, 
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Note that the overall factor of 1/2 for creating a pair of 
identical particles is canceled by P ilv P^ v = QuvQ^ v = 2. 



IV. BENCHMARK COLOR OCTET MODEL 

We now specialize our results to our benchmark model: 
a single real, color-octet, electroweak neutral scalar 

s a us EM. 

C S = \{D,S a ) 2 - \m 2 (S a ) 2 - *(S a ) 2 H^H 
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Here we have included one Higgs-portal interaction, with 
coupling well as additional renormalizable self- 

interactions among the S a - a quartic interaction with 
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coupling u), and a cubic interaction with (dimensionful) 
coupling us. The cubic interaction implies the S a can 
decay into two gluons through a triangle loop. Since the 
Lagrangian possesses a Z2 symmetry when /i s = 0, small 
values of fis a re technically natural. This allows us to 
assume that fis is large enough to cause prompt S a de- 
cays, but is otherwise small enough to not play a role in 
our calculations of di-Higgs production. 

The effect on single Higgs production of light scalars 
through the Higgs-portal interaction was precisely the 
subject of our Ref. [12 . There we showed that single 
Higgs production could be dramatically reduced for a 
wide range of negative K, and scalar masses in the range 
Ms ~ 100-300 GcV. The reduction comes from de- 
structive interference between the top quark loop and 
the loop of S a . The suppression could reach 75 — 90% 
without excessive tuning. For the most extreme cancel- 
lation, relatively large couplings k ~ — 1 and light scalars 
Ms < 200 GeV were required. Taking even lighter scalar 
masses, the scalar loop contributions can overwhelm the 
top quark contribution, and the single Higgs production 
cross section can again be 0(1) times the Standard Model 
value. 

As was discussed in Ref. [12], color octet scalars are 
dominantly pair-produced and decay (via a loop of S a ) 
to a pair of gluons. Thus, the signature is two pairs of jets 
with equal invariant mass, making four total jets. This 
signature provides a handle to help distinguish it from 
the multi-jet background. However, when the scalars 
are light, they create relatively soft jets that are not ef- 
ficiently triggered on due to the large background. In 
addition to the large rate and combinatorial hurdles, ex- 
tracting the signal is further complicated by pileup. A 
comparison between ATLAS [37] and CMS [3SJ multi- 
jet searches demonstrates the difficulties in placing con- 
straints in higher luminosity runs. The ATLAS [37] 
search was able to disfavor real color octet scalar between 
100-125 GeV. The low reach in Ms was possible because 
this study was performed with 36 pb -1 of data where 
the luminosity was sufficiently low that trigger thresholds 
could be kept small. The CMS 35J search, by contrast, 
was performed on 1 fb _1 of data, where the jet trigger 
threshold was 70 GeV. This prevented the CMS study 
to be sensitive to objects lighter than 300 GeVQ By 
adopting specialized pre-scaled triggers at lower jet pr, 
improvements to the low mass octet search may be pos- 
sible, but to the best of our knowledge, no such studies 
have been completed. 

It is also important to consider the possible range of the 
dimensionless Higgs-portal coupling k, particularly for 
negative values. Large negative k could destabilize the 
scalar potential at large field values. The naive tree-level 
bound can be obtained by rewriting the quartic interac- 
tions of the scalar potential as (y/u}(S a ) 2 — y/XhWH) 2 + 



The jet px cut was 150 GeV in this study. 



2y/w\i t S%Hi H , and thus | k\ < [lOj where A/i is 

the Higgs quartic coupling. However, in Appendix |C] we 
calculate the renormalization group equations for k, w, 
and Xh, showing that negative K runs rapidly to smaller 
(absolute) values at a larger renormalization scale. This 
suggests that the bounds on negative k are expected to be 
significantly relaxed once one uses the renormalization- 
group improved potential. 



V. ENHANCED DI-HIGGS PRODUCTION 

We now calculate the di-Higgs production cross section 
for two distinct scenarios: mh = 125 GeV, and nih > 
125 GeV. In the case rrih = 125 GeV, we consider both 
negative and positive k, since a wide range of (k,Ms) 
parameter space is allowed by a (generous) interpretation 
of the 125 GeV resonance observed by LHC as a Higgs 
boson. For rrih > 125 GeV, we consider only negative k, 
since this is the only region of (k, M s ) space where the 
(heavy) single Higgs production can be suppressed below 
the LHC bounds [H]. 

A. m h = 125 GeV 

We are finally in a position to evaluate the di-Higgs 
production cross section, Eq. ( |18[ ), comparing to the 
Standard Model value while scanning over the scalar pa- 
rameters Ms and k. In Fig. [3] we show numerical results 
for the ratio a(pp -> /i/i)top+scaiar/c(pp -> hh) top as a 
function of M s and K for LHC at = 8 TeV. The 
scalar and top calculations were performed at leading or- 
der; however by taking a ratio to the Standard Model 
result, higher order corrections in a s to di-Higgs produc- 
tion should largely cancel out. Fig.|3]contains three plots: 
one with positive and negative k with Ms < 400 GeV; 
and two "zoomed in" plots that more easily demonstrate 
the large effects at smaller Ms for both positive and neg- 
ative k. These contours are overlaid onto colored re- 
gions that show the single Higgs production ratio cr(pp — ¥ 
fr)top+scaiar/o-(pp -> ft) top j derived from Ref. [12]. We see 
that once Ms ^ rrih, the di-Higgs production contours 
follow the single Higgs production region shapes fairly 
well. This was partly anticipated from the coefficients of 
the effective operators, Eqs. (|7|)>P|> that were found to be 
the same for h/vG^GW and tY/(2v 2 )G° v G a " u . Inter- 
estingly, the magnitude of the di-Higgs production cross 
section is not estimated particularly well from the effec- 
tive operators. In Appendix [B] we compare the results 
for the exact momentum-dependent one-loop calculation 
and the effective effective operators. This we do by show- 
ing the negative k "zoomed in" plot identical to Fig. |3jc) 
overlaid with the effective operator results. It is perhaps 
not surprising that these results do not agree, since the 
effective operators neglect momentum-dependent correc- 
tions of order s/M 2 that is order one or larger in the re- 
gion we show. Perhaps more interesting, however, is that 
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the effective operators do reproduce the shape of the con- 
tours of the di-Higgs cross section ratio in (ft, Ms) when 
the ratio a(pp -> hh) sca x aT+top /a(pp -> hh) top itself is 
large (say, > 10) and Ms ^ rrih- Of course in the region 
Ms < rrih, as well as when the ratio itself is relatively 
small (say, < 10), the effective operators do not repro- 
duce either the magnitude or shape of the cross section 
ratio results. 

If we generously allow the recent LHC results to sug- 
gest 0.5 < a(pp -> /l) SC alar+top/o'(pp -> /l)top < 2, then 

there are three regions of interest: 

• Negative re, below and to the left of the red re- 
gion. This region has the largest di-Higgs cross 
section enhancement, since it is the region with 
largest negative re for a given Ms- In this region 
the scalar contributions to the single Higgs produc- 
tion are roughly twice the size, but opposite in sign, 
to the top contribution, resulting in a single Higgs 
production rate that can be smaller, the same as, 
or larger than the SM single Higgs production rate. 
The di-Higgs cross section, by contrast, is enormous 
- between several roughly 100-1000 times the Stan- 
dard Model cross section. 

• Negative re, above and to the right of the red region. 
This region has a single Higgs production cross sec- 
tion that is smaller than the Standard Model. Nev- 
ertheless, the di-Higgs production cross section can 
be up to about 10 times the Standard Model rate. 
In this region there is a strong correlation between 
(slight) suppression of the single Higgs production 
rate and the (much larger) enhancement of the di- 
Higgs production rate. 

• Positive re. In this region, the single Higgs produc- 
tion rate is larger than the Standard Model. Re- 
stricting to a(pp -t h) sca i ar+ t op /o-(pp -t h) top < 2, 
we see the di-Higgs enhancement is negligible. It is 
in this region that smaller positive re leads to con- 
structive interference for single Higgs production 
and destructive interference for di-Higgs produc- 
tion. 



B. m h > 125 GeV 

Consider now the case where the 125 GeV particle ob- 
served by LHC is, in fact, an imposter. As we showed 
in Ref. |12j . there is viable parameter region where the 
Higgs can be much heavier, and yet, be safe against the 
LHC search bounds on single Higgs production. This oc- 
curs when the colored scalar contributions to single Higgs 
production interfere destructively with the top contribu- 
tion, leading to a large region of highly suppressed single 
Higgs production rate at LHC. We show this in Fig. [I] 
where the viable region lies entirely in the red and some 
of the blue region, where the single Higgs cross section 
at LHC, a(pp -> h) sca _i ar+top /a(pp -t h) top < 0.1-0.25. 



Here we see that there is a perfect correlation between 
suppression of single Higgs production and enhancement 
of di-Higgs production. This occurs for all of the Higgs 
masses shown, with the largest di-Higgs enhancement oc- 
curring for the smaller Higgs masses, rrih = 160, 200 GeV, 
large negative re, and Ms 5s 250 GeV. Here the enhance- 
ment is between about 5 to 200 times the Standard Model 
rate (for the given Higgs mass) . This is the smoking gun 
for a heavy hidden Higgs boson in the case where the 
125 GeV particle is an imposter. 



C. Kinematic Distributions of Enhanced Di-Higgs 
Production 

In addition the large increase in the di-Higgs produc- 
tion rate, there are also modifications to the kinematical 
distributions of di-Higgs production. Here we turn to 
Monte Carlo to simulate the signal. We implemented 
the gg —> hh processes into MadGrapM [5U] by mod- 
ifying the necessary HELAS [5TJ [52] routines, including 
both Lorentz structures (Puv & n d Q^v) and retained the 
full momentum dependence. To evaluate the Passarino- 
Veltman one-loop functions, the gg — » h, gg — > hh ampli- 
tudes were interfaced with the LoopTools [53] package. 
We utilize CTEQ6L1 parton distribution functions with 
the scale choice /if — M-R = 2 x for all simulations. 

We first look to basic kinematic distributions of the 
Higgs bosons. The pr and rapidity (Y^) spectra are 
shown below in Fig. [5] for several choices of Ms; the 
distributions are area-normalized to focus on the shape 
difference. Note that the PT,h spectrum peaks roughly at 
the mass of the particle dominating the loops, near rrit in 
the SM and Ms in all non-SM cases. For the SM case this 
feature can be traced to an enhancement in the diagrams 
when s ~ ml [T]. A similar enhancement occurs in the 
scalar loops at s ~ Mj. However, unlike in the SM where 
the triangle and box diagrams always destructively inter- 
fere, the interference among the various scalar diagrams 
depends on the sign of re, as we have seen. For positive re, 
the scalar boxes and triangles also interfere destructively. 
The severity of the interference depends on |re|, since the 
box contributions grow as 0(re 2 ), whereas the triangles, 
O(re). 

In Fig. [5] the Higgs-portal coupling was fixed to re = 
— 1, while the scalar mass was varied. However, as there 
are both O(re) and 0(k 2 ) contributions to gg — > hh (at 
amplitude level), each weighted by different kinematic 
functions, the pr,h spectrum does carry some re depen- 
dence. As we can see from Fig. [6] which shows the over- 
layed PT,h spectra for three different re values and fixed 
Ms — 150 GeV, the re dependence is fairly small. 

Next, we come to the differences between using effec- 
tive operators and retaining the full one-loop momentum- 
dependence. In Fig. [7J we take Ms = 150 GeV, re = — 1 
and compare the pt,k distributions between using the ef- 
fective operators [Eqs. (p~3|)- ( 14)] against the one-loop re- 
sults [Eqs. tfTBft-dlTb 



Clearly, the distributions are qual- 
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FIG. 3. Di-Higgs as well as single Higgs production cross sections normalized to the Standard Model values. The solid lines 
show a(pp — > ft/i) SC aiar+top/c r (pp — >• hh) top , while the colored regions show when a(pp — > h) Bca \ aT+top /a(pp — > h) top is less than 
the numerical value labeling each region. Panels (b) and (c) are "zoomed in" versions of panel (a), to more clearly see the 
contours in the light Ms < 200 GeV region. 



itatively distinct throughout the pT,h range. For the ef- 
fective operator calculation, the high-p^ tail is governed 
by kinematics alone, where as there is an extra suppres- 
sion from the form factors once we incorporate momen- 
tum dependence. An accurate px spectra is vital for 
phenomenology; the Higgs px will be transferred to its 



decay remnants, and the details of the remnant kinemat- 
ics are a necessary ingredient for successfully separating 
signal from SM background, regardless of the identity of 
the final state particles. 
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FIG. 4. Same as Fig. [3] but for heavier Higgs masses (assuming the 125 GeV particle observed by LHC is an imposter). The 
viable parameter space is roughly in the red and blue regions where a(pp — > /i)scaiar+to P /o"(w ~> Mtop Si 0.25. Whenever 
rrih > 2 Ms,2mt there is little space for a suppressed o(jpp — > h). This occurs because the loop contributions to the amplitude 
become complex, and cancellation must be arranged simultaneously between the real and imaginary parts of the top quark and 
scalar loops. 



VI. DISCOVERING DI-HIGGS 



The signals and discovery strategies for di-Higgs pro- 
duction crucially depend on the mass of the Higgs boson 
and how much the rate is enhanced. We will not at- 
tempt to cover all scenarios, and instead focus on two 
distinct cases: = 125 GeV and = 200 GeV. 

The rrih — 200 GeV scenario is fairly representative of 



a generic heavier Higgs mass (m/j > 200 GeV), since the 
branching fractions into massive gauge bosons dominate. 
We will sketch the best signals for each case and consider 
strategies to improve the discovery potential. A detailed 
study of the full background rates and the optimal cut 
strategy for a given Higgs mass is beyond the scope of 
this work. 

In the case — 125 GeV (where the particle ob- 
served by LHC is taken to be the Higgs boson), the rate 
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(a) (b) 




FIG. 5. Kinematic distributions for the Pt,h of the Higgs boson in di-Higgs production at LHC. Panel (a) shows (l/a)da(gg — > 
hh)/dpT.h for the SM alone (red) and the SM combined with loops of scalars, Ms = 160 GeV (green), Ms = 200 GeV (orange), 
Ms = 250 GeV (gray) and Ms = 300 GeV (cyan). The Higgs-portal has been taken to be re = — 1 for all cases. The right panel 
shows the Higgs rapidity distribution (l/a)da(gg — > hh)/dYh for the same scalar values. 




FIG. 6. Comparison between the Pt,h spectra in gg — > hh 
for three different values of re. The scalar mass is held fixed 
at 160 GeV but we vary the Higgs-portal coupling; k = —1.0 
(green), k = —1.5 (blue), re = —0.5 (magenta). The distribu- 
tions are area normalized to focus on the shape differences. 



for single Higgs production is constrained to be roughly 
the Standard Model value. As we saw from Sec. |V A| 
given a range in the single Higgs production rate, there 



is a corresponding range of re, Ms parameter space. The 
wider the range, the larger the variation in the di-Higgs 
production rate. 

Here we are interested in the absolute rate of di-Higgs 
production. While we have not calculated beyond one- 
loop, Ref. [3] found the K factor for di-Higgs production 
in the Standard Model to be 0(2). This means <r(pp — > 
h)/a(pp -> hh) ~ 2500 at y/s = 8 TeV. Because there 
are two Higgs bosons in the final state, the rate for di- 
Higgs to feed into any characteristic single-Higgs final 
state (77 + X, ZZ* + X, etc.) is obviously doubled. We 
expect that these di-Higgs events would be captured by 
the inclusive single Higgs searches at the LHC. In the 
absence of an observation, this should allow ATLAS and 
CMS to obtain an estimate of the upper bound on the di- 
Higgs enhancement. Dedicated analyses of existing y/s — 
7 and 8 TeV data that focus on the exact di-Higgs final 
state and kinematics, such as 77 + 66, would undoubtedly 
increase the sensitivity to the signal. 

The discovery prospects for some di-Higgs final states, 
including with enhanced production rate, have been out- 
lined in Ref. [2TM21 US [HI [35] . However, in the time 
since those studies there have been significant advances 
in the usage of jet substructure as a discovery tool, for 
both hadronic (66) [53H55] and r + r~ [571 [5B] resonances. 
Given the relatively high px that the Higgs bosons can 
carry, i.e., Fig. [5] we expect di-Higgs events to be well 
suited to substructure techniques. 

The rates, in fb, for several potential di-Higgs final 
states are shown in Table |U All rates are functions of the 




FIG. 7. The difference in the pr,h spectrum of the Higgs 
bosons from di-Higgs production between using the full one- 
loop momentum dependence (solid line) versus the effective 
operators (dashed line). In this figure we have taken Ms = 
150 GeV, k — — 1 . The distributions are area normalized to 
focus on the shape differences. 





final state 


rates in fb 
8 TeV 14 TeV 


m h = 125 GeV 


77 + bb 
t + t~ + bb 


0.019V 0.084V 
0.53 X 2.4 X 
0.029V 0.13V 


m h = 200 GeV 


w+ww+w- 
zzzz 
w+w-zz 


1.21X 5.73V 
O.UX 0.68V 
0.42X 1.97X 



TABLE I. The di-Higgs signal production rates into the dif- 
ferent modes when mi = 125 GeV and rrih = 200 GeV. All 
rates, which include a K factor of 2.0, must be multiplied 
by the enhancement factor X which can be read off from 
the di-Higgs contours in Fig. [3] [rrih = 125 GeV) and Fig. [4] 
(m h = 200 GeV). 



enhancement X and we have taken the JGfactor to be 2.0 
throughout. We find that varying the parton distribution 
set changes the cross-section by O(20%). 

Now, in the alternative case rrih = 200 GeV, the single 
Higgs production cross section is constrained to be less 
than the present LHC bounds. This basically requires 
us to choose (k,Ms) within the red or blue region of 
Fig. [4] By taking rrih — 200 GeV, both diboson decay 
modes are open, i.e., rrih > 2m^, 2rriz, and the di-Higgs 
production phenomenology becomes fairly insensitive to 
the exact value of the Higgs mass. In Table [I] we show 
the rates for the various combinations of dibosons. The 
hh —> W + W~W + W~ is the largest, and provides several 
opportunities involving same-sign dileptons or trileptons 



plus missing energy. Here again, jet substructure tech- 
niques could provide valuable additional sensitivity given 
the large px of the Higgs bosons in di-Higgs production. 



VII. DISCUSSION 

We have demonstrated that di-Higgs production at the 
LHC could be many orders of magnitude larger than the 
Standard Model in a the presence of light colored scalars. 
The large enhancements are possible if rrih = 125 GeV, 
consistent with the recent data from LHC, as well as if 
rrih > 125 GeV, should the 125 GeV particle turn out 
to be a Higgs imposter. The latter region is natural 
in a model with light colored scalars, since this model 
was already demonstrated to effectively hide the Higgs 
below the LHC limits [12] . The largest effects occur 
when the Higgs-portal coupling re is negative and col- 
ored scalars that are lighter than the electroweak break- 
ing scale. Our detailed numerical results were performed 
in a model with a single, real, color octet scalar. Di- 
rect production of the colored scalar and decay to pairs 
of gluons is presently constrained by LHC data only for 
Ms 5s 125 GeV. Hence, a wide range of (re, Ms) space 
remains viable that can lead to the huge di-Higgs en- 
hancement rates that we find. 

Our calculation of the production rates were performed 
for a general set of colored scalars at leading order, and 
thus up to the overall multiplicity, are independent of 
the electroweak quantum numbers. The detailed branch- 
ing fractions of the Higgs will, of course, depend on the 
electroweak quantum numbers. Generally this leads to 
O(l) effects on the branching fractions, whereas the di- 
Higgs production rate enhancements we found can be 
orders of magnitude larger. It would be interesting to 
perform a more detailed investigation of the correla- 
tion between single Higgs production and decay, with 
other representations of colored scalars, versus the di- 
Higgs production cross section [TTJ [T31 [551 ISO] ■ In ad- 
dition, electroweak charged scalars will alter a wider set 
of processes than electroweak singlets: pp — > VV, where 
V = 7, Z°, W will all change due to loops of electroweak 
charged scalars, in addition to pp — > hh. Restricting to 
electroweak neutral colored scalars, we were somewhat 
forced to use octets. This is because smaller representa- 
tions, for instance the color triplet or sextet, are stable 
when the scalars are electroweak singlets. 

The signals of di-Higgs production provide several ex- 
citing opportunities for the LHC. We have already 
seen that di-Higgs production could have already con- 
tributed a fraction of the single Higgs production rate, 
with the second Higgs missed (or unidentified) in the 
events. If rrih = 125 GeV, among the more interest- 
ing di-Higgs signals includes hh — > 7766, hh — > T + r~bb, 
and hh — > t + t~t + t~ . Each of these signals has two 
pairs of particles that reconstruct to two Higgs bosons. 
If rrih > 125 GeV, there are several channels resulting 
from hh — > W + W~W + W~ , including same-sign leptons 
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plus missing energy as well as trileptons plus missing en- 
ergy that could be seen at LHC. The latter two searches 
might first obtain evidence from recasted supersymmetric 
searches, while a more dedicated search strategy would 
undoubtedly improve the sensitivity. Jet substructure 
may also provide a valuable additional strategy for both 
low and high Higgs masses. 



Appendix A: Passarino-Veltman One-Loop 
Functions 



The Passarino-Veltman functions are defined by: 



C Q (pi,pj : m) = 
Do(pi,p 3 ,Pk ■ m) = 



in {q 2 — m 2 )((q + pi) 2 — m 2 )((q + Pi + Pj) 2 — m 2 ) 
d 4 q 1 



in {q 2 - m 2 )((q + p^ 2 - m 2 )((q+p l + pj) 2 - m 2 )((q + pi + pj + p k ) 2 - m 2 ) 



(Al) 
(A2) 



Appendix B: Effective Couplings versus One-loop 
Results 



Comparison of exact one-loop result against that ob- 
tained from using the effective couplings in Fig. [8] Here 
we only show a "zoomed in" region for small Ms and 
negative k, where the differences between the calcula- 
tions are largest. Nevertheless, we have verified that for 
Mi 3> v,mh, the contours asymptotically agree. 
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FIG. 8. Comparison of exact one-loop result against that 
obtained from using the effective couplings 



Appendix C: Renormalization Group Equations 

The renormalization group equations for k, uj, and 
(taking [is to be small and negligible) are: 



16n 2 



= 2k 2 -9Kg 2 c +6K\ h + 10 kuj (CI) 



16tt 2 ^ = k 2 + 24w 2 - l%g 2 c uj + 12 g* (C2) 

where g c is the 5*17(3) coupling and Xh is the Standard 
Model Higgs quartic coupling. All terms in the n renor- 
malization group equation are proportional to k, since 
k = is technically natural. The running of Xh is also 
changed by the Higgs portal interaction: 



16tH 



, dXh 
lit 



(C3) 



where j3\ h is the /3-function for the Higgs quartic coupling 
coming from SM interactions. 



For large positive k, Eq. (CI) shows that n will grow 



and ultimately hit a Landau pole, similar to what hap- 
pens with the H iggs quartic coupling. For large negative 
k, however, Eq. ( CI ) shows that the Higgs-portal interac- 



tion rapidly decreases until the other terms in the renor- 
malization group equation become important. Hence, 
a one-loop renormalization group improved potential is 
necessary to analyze the stability of the effective poten- 
tial in the presence of a large negative n coupling. Our 
estimates suggest n > —2 is perfectly acceptable, how- 
ever a detailed analysis of the scalar potential is beyond 
the scope of this paper. 
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